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ABSTRACT 

We report the discovery of 112-ms X-ray pulsations from RX J0822— 4300, the compact central 
object (CCO) in the supernova remnant Puppis A, in two archival Newton X-Ray Multi-Mirror Mission 
observations taken in 2001. The sinusoidal light curve has a pulsed fraction of 11% with an abrupt 
180° change in phase at 1.2 keV. The observed phase shift and modulation are likely the result of 
emission from opposing thermal hot spots of distinct temperatures. Phase-resolved spectra reveal an 
emission feature at En ne = 0.8 keV associated with the cooler region, possibly due to an electron 
cyclotron resonance effect similar to that seen in the spectrum of the CCO pulsar IE 1207.4—5209. 
No change in the spin period of PSR J0821— 4300 is detected in 7 months, with a 2a upper limit on 
the period derivative of P < 8.3 x 10~ 15 . This implies limits on the spin-down energy loss rate of 
E < 2.3 x 10 35 erg s _1 , the surface magnetic dipole field strength B s < 9.8 x 10 11 G, and the spin- 
down age r c > 220 kyr. The latter is much longer than the SNR age, indicating that PSR J0821— 4300 
was born spinning near its present period. Its properties are remarkably similar to those of the two 
other known CCO pulsars, demonstrating the existence of a class of neutron stars born with weak 
magnetic fields related to a slow original spin. These results are also of importance in understanding 
the extreme transverse velocity of PSR J0821— 4300, favoring the hydrodynamic instability mechanism 
in the supernova explosion. 

Subject headings: stars: neutron — pulsars: (RX J0822-4300, PSR J0821-4300, IE 1207.4-5209, 
PSR J1210-5226, PSR J1852+0040) — supernova remnants: individual (Pup- 
pis A) 
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1. INTRODUCTION 

Recent timing studies of central compact objects 
(CCOs) in supernova remnants (SNRs) have found young 
neutron stars (NSs) whose spin-do wn is imperceptible 
due t o a weak dipole magnetic field (|Gotthelf fc Halpernl 
|2008|) . More generally, CCOs are defined by their steady, 
predominantly thermal X-ray emission, lack of optical or 
radio counterpa rts, and absence o f puls ar wind nebula e 
(see reviews by iPavlov et all 120041 and iDe Lucal [20081 ) . 
Of the six firm members of the CCO class, two are 
identified explicitly as pulsars, while searches for pulsa- 
tions from the other CCOs have been unsuccessful. In 
particular, claims have been made for several candidate 
spin periods of RX J0822-4300, the CC O in Puppis A 
([Pavlov et alJl!999t iHui fc Beckerll2006an . but none has 
been confirmed. 

RX J0822— 4300 has long been considered the stellar 
remnant of the supernova explosion that produced Pup- 
pis A, given its central location, lack of an o ptical coun- 
terpart, and hot thermal spectral properties (jPetre et al.l 
119821 1 1 996( 1. Puppis A itself is an oxygen-rich SNR, 
similar to Cas A. Its age is estimated as 3.7 ± 0.4 kyr 
from the proper motion of its oxygen knots (Winkler 
et al. 1998), and its distance d = 2.2 ± 0.3 kpc is de - 
rived from H I absorption features (jReynoso et al.i r2002) . 
More recently, an unusually large transverse velocity 
of sa 1600 km s~* was measured for RX J08 22-4300 
(jWinkler fc Petrel [2007t IHui fc Becked l2006bh . testing 
proposed mechanisms of natal kick. 

In this Letter we present a reanalysis of a pair of timing 
observation of RX J0822— 4300 acquired with the New- 



ton X-Ray Multi-Mirror Mission (XMM-Newton) ob- 
servatory. We find a highly significant pulsed signal 
in both data sets that imposes an interesting upper 
limit on the period derivative. PSR J0821-4300 be- 
comes the third CCO pulsar, with timing and spectral 
properties similar to those of the t wo other CCO pul- 
sars, PSR J1852+0040 in Kes 79 (|Gotthelf fc Halpernl 
[20051: lHalnern et all [2007h and PSR J121 0-5226 in 
PKS 1209-51/52 (jGotthelf fc Halperd[200l . We dis- 
cuss the physical basis for the class of CCO pulsars and 
how their natal properties differ from radio pulsars and 
magnetars. 

2. XMM-NEWTON OBSERVATIONS 

Two archival XMM-Newton observations of RX 
J0822— 4300 are best suited to search for coherent pul- 
sations. These were obtained on 2001 April 15 and 
Novem ber 8 using the Euro pean Photon Imaging Camera 
(EPIC: [Turner et al.ll2003T) operating in "small- window" 
mode (4f3x4f3 field-of-view) . This mode provides 5.7 ms 
time resolution, allowing a search for even the most 
rapidly rotating young pulsar. The EPIC pn detector is 
sensitive to X-rays in the nominal 0.1 — 12 keV range 
with energy resolution AE/E w 0.1/-\/£ , (keV). The 
medium filter was used with the target placed at the 
default EPIC pn focal plane location for a point source. 
The concurrent EPIC MOS data are not used in this 
work because of its more limited (2.7 s) time resolution. 

We reprocessed the EPIC pn data using the Science 
Analysis System version SAS 7.0.0 (2006028.1801) and 
screened the photon event lists using the standard cri- 
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Fig. 1.— Discovery of PSR J0821-4300 in Puppis A using XMM-Newton EPIC pn data acquired on 2001 April 15 (left) and November 
8 (right). The power spectrum of 1.5 — 4.5 keV photons is extracted from a 30" radius aperture at the location of RX J0822— 4300. The 
detection threshold for a blind search for a period P > 12 ms is indicated. Inset: Background subtracted pulse profiles of PSR J0821— 4300 
in two energy bands illustrating the phase shift between the two; the soft photons are extracted from a 18" radius aperture to minimize 
the SNR background contamination. 



TABLE 1 

XMM-Newton Timing Results for PSR J0821-4300 



Epoch Dur. Bkgnd. a Source b Period c f v d 

(MJD) (s) (s- 1 ) (s- 1 ) (s) (%) 



52,014.325 22,637 0.131(4) 0.653(8) 0.112799416(51) 11(1) 
52,221.765 22,511 0.158(5) 0.668(8) 0.112799437(40) 11(1) 



(|Buccheri et al.lll983l ) and localized the pulsed emission 
with a peak statistic of Z\ = 48.94 and Z\ = 50.20 for 
the 2001 April and November observations, respectively. 
These statistics correspond to 99.991% and 99.995% con- 
fidence, respectively, after allowing for the number of in- 
dependent trials (2AT span /P mm ) in a blind search for 
-fmin > 12 ms, the Nyquist limit. The timing results are 
summarized in Table [U 

The periods derived from the Z\ test are statistically 
identical. The periodogram derived from each data set is 
shown in Figure Q] along with the pulse profile folded at 
the best period. With a single coherent fold using both 
data sets, a total power of Z\ = 95.98 is recovered with- 
out the need for a period derivative. This corresponds to 
a negligible probability of chance occurrence. The formal 
result is P — (1.2±3.6) x 10~ 15 s s — x . The uncertainty is 
lcr and is computed by propagating the uncertainties on 
the individual period measurements. From the 2a upper- 
limit of P < 8.3 x 10~ 15 s s _1 we constrain the spin- 
down power of PSR J0821-4300 to E = 4ir 2 IP/P 3 < 
2.3 x 10 35 erg s _1 , the surface dipole magnetic field 

strength to B s = 3.2 x 10 19 VPP < 9.8 x 10 11 G, and 
the characteristic age to r c = P/2P > 220 kyr. 

The background-subtracted pulsed fraction is 11 ± 1% 
for both observations. The pulsed fraction is defined here 
as f p = j'V(pulsed) /N (total), where we choose the min- 
imum of the folded light curves as the unpulsed level. 
The pulse shape is evidently unchanged between obser- 
vations, suggesting a stable underlying emission process. 
No other significant signal is detected in either data set. 
In particular, we examined the two we aker candidate sig- 
nals reported bv lHui fc Beckerl (|2006aD from these obser- 
vations. Although we see these peaks, neither is statis- 
tically significant or repeatable between the two obser- 
vations. We also analyzed all archival X-ray data sets 
that might yield detections of PSR J0821-4300 and fur- 
ther constrain its timing parameters, but none proved 
sensitive enough. 

A search for the signal in the softer energy band re- 
veals that the pulse is indeed present, below 1.0 keV, 
but shifted in phase by half a cycle relative to the hard 



Note. — Statistical uncertainty in the last digit is given in 
parentheses. 

a Background rate in the 1.5 — 4.5 keV band in a 0'.5 < r < 0'8 
annular aperture, corrected for dead-time. b Count rate in the 
1.5 — 4.5 keV band in a r = 0f5 aperture, corrected for dead- 
time and background. Period derived from a Zj test. Uncer- 
tainty is lcr, computed by the Monte Carlo method described 
in IGotth elf ct al. (1999). d Observed pulsed fraction, defined as 
f p = A r (pulsed)/A r (total), after subtracting background. 

teria. Both observations were uncontaminated by flare 
events, providing w 22.5 ks of good EPIC pn expo- 
sure time during each epoch. After taking into ac- 
count the CCD readout dead-time (29%), this trans- 
lates to 16 ks of live-time per observation in "small 
window" mode. Photon arrival times were converted 
to the solar system barycenter using source coordinates 
08 h 21 m 57;37, -43°00'17'/0 (J2000.0) derived from a con- 
temporaneous Chandra observation. In the following 
analysis we choose a source aperture 0'5 in radius cen- 
tered on RX J0822— 4300 to maximize signal over back- 
ground, except where otherwise noted. This region en- 
closes > 80% of the source flux. Count rates and expo- 
sure times are presented in Table [TJ 

2.1. Timing Analysis 

To search for a pulsed signal, arrival times of photons 
were initially extracted in the energy band 1.5 — 4.5 keV 
under the hypothesis that hard X-rays coming from a 
smaller area might be more strongly pulsed than soft X- 
rays from the full stellar surface. A 2 23 -bin Fast Fourier 
Transform was used. The most significant signal de- 
tected was P = 112 ms in both EPIC pn data sets, with 
no higher harmonics. We constructed a periodogram 
centered on this signal using the Z\ (Rayleigh) test 
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Fig. 2. — Pulse phase as a function of energy for PSR 
J0821— 4300. The phases and errors are determined by cross- 
correlating energy-selected folded light curves with a master tem- 
plate. The phase is computed for overlapping 0.5 keV wide en- 
ergy slices incremented in 0.1 keV steps. Phase zero is defined as 
the minimum of the modulation as shown in Figure [T] The re- 
sult using the April 15 data set is shown; the November 8 data 
yields an equivalent result. Both show a half cycle phase shift at 
Eshift = 1.1 - 1.2 keV. 

pulse. Figure [2] shows that the phase shift is remarkably 
abrupt, with a half cycle step (A0 = 0.5) at an energy of 
1.1 — 1.2 keV. The location of the shift is unresolved to 
better than ~ 0.25 keV, limited by the counting statistics 
and the EPIC pn energy resolution of ~ 0.1 keV at 1 keV. 
This shift effectively cancels the signal in a full-spectrum 
search. 

2.2. Phase- Resolved Spectroscopy 

iHui fe Becker! (|2006af ) performed a comprehensive 
spectral analysis of the XMM-Newton data on RX 
J0822— 4300, showing that a two-component fit is nec- 
essary. They preferred a double blackbody model, which 
is also compatible with our expectation that surface ther- 
mal emission dominates over any other source of X-rays 
from CCOs. We can investigate the cause of the phase 
shift in the pulse profiles by examining pulse phase- 
resolved spectra and searching for features that may be 
phase-dependent. In particular, we extracted two spec- 
tra, one in the phase band centered on the pulse peak 
found in the harder energy band (<f> = 0.5), as shown 
in Figure [1] and the other centered on the phase band 
found for the softer emission (</> = 0.0). Each spectrum 
covers a phase interval 0.4 cycles and is referred to herein 
as the hard and soft phase spectra. 

We combined data from both observations and grouped 
the two phase resolved spectra with a minimum of 100 
counts per channel. These were fitted to the double 
blackbody model in the 0.4 — 4.0 keV energy range. We 
find that both spectra can be fitted with the same com- 
ponent temperatures, with phase-dependent normaliza- 
tion. The fitted spectra are shown in Figure [3] and the 
best fit parameters are listed in Table [U For the soft 
phase there is clear evidence of positive residuals from 
the continuum model in the 0.7 — 0.9 keV range. This 
is well modeled by the addition of a Gaussian emis- 
sion line of -Eiinc = 0.79 ± 0.02 keV with an equivalent 
width of W ~ 55 eV. Comparing a fit without and with 
the line feature results in a change from x 2 — 69.41 
to x 2 = 53.45, for 70 and 67 degrees-of-freedom, re- 
spectively. For the hard-phase spectrum, the continuum 
model provides a satisfactory fit without the need for an 
added line component (W <J 9 eV). 

Given the excellent fits to the double blackbody model 



TABLE 2 

Phase Resolved XMM-Newton Spectral Results for PSR 
J0821-4300 



Parameter 


Soft Phase 


Hard Phase 


N H (cm" 2 ) 
feTi (keV) 
kT 2 (keV) 
L^BB) (ergs- 1 ) 
L 2 {BB) (ergs" 1 ) 
Ai(BB) (cm 2 ) 
A 2 (BB) (cm 2 ) 


4.8 ±0.06 

o.2i+8:°g 

0.42l«;° 4 
4.1 x 10 33 
2.5 x 10 33 
2.1 x 10 12 
7.8 x 10 10 


(Linked) 
(Linked) 
(Linked) 
3.3 X 10 33 
3.2 x 10 33 
1.7 x 10 12 
9.9 x 10 10 


Gaussian Line Component 


E line (keV) 
"line (eV) 
W (eV) 


0.79 ± 0.03 

cq+33 
JO -22 

« 55 


(Linked) 
(Linked) 
< 9 



X z (dof) 125.6 (153) 



Note. — Spectral fits in the 0.4 — 4 keV energy band, for two 
phase regions, 0.4 cycles wide, one centered on the pulse peak found 
in the softer energy band ( "soft phase" ) and the other centered on 
phase found in the harder energy band ("hard phase"). All errors 
are at the 90% confidence level. 

the most natural explanation of the emission from PSR 
J082 1—4300 is that of two discrete "hot spots" on the 
surface of the neutron star, corresponding to the two 
blackbody temperatures and areas. As the star rotates, 
the hard emission dominates over the soft component 
when the hotter side is facing the viewer, while half a 
cycle later, the soft emission dominates as the harder 
emission is now (partially) eclipsed by the star. Based 
on their phase relationship, evidently 180°, the two spots 
are on opposite sides of the NS. 

Several measurements favor this explanation: 1) both 
the soft and hard phase spectra can be fitted to a dou- 
ble blackbody model with the same component tempera- 
tures, but phase dependent normalizations, 2) the phase 
shift energy (-E s mft = 1.1 — 1.2 keV) is in agreement with 
the cross-over energy between blackbody components for 
the above best fit model (E cmss = 1.16 keV; cf. Figure [2] 
and Figure[3]), 3) for this model fit, the magnitude of the 
pulse modulation is consistent with the change of the de- 
rived areas for both the soft and hard phase spectra (see 
Table EJ). 

In this picture, the geometry is strongly constrained. 
The small modulation requires that the spin axis be 
nearly aligned with the line-of-sight. Furthermore, the 
hot spot axis must be nearly perpendicular to the spin 
axis since both temperature components are evident 
throughout the rotation phase. Future modeling of the 
observed spectrum and light curves will allow us to 
test this hypothesis, of opposed thermal hot-spots for 
PSR J0821— 4300, and quantify the hot-sp ot and view- 
ing geometries (e.g., iPerna fc GoTt helf 2008h . This study 
should also provide important constraints on the pre vi- 
ously allowed spectral models ()Hui fc Beckerl l2006al) of 
the pulsar. 

3. DISCUSSION 

3.1. Physical Basis of the CCO Class 

In the dipole spin-down formalism, the 2a lower limit 
of t c > 220 kyr on the characteristic age of PSR 
J082 1—4300, in combination with its true (SNR) age of 
T = 3.7 kyr, implies that its initial spin period Pq = 



4 



Gotthelf & Halpern 



- PSR J082 1-4300 



2BB+GAU FIT 



= 0.0 




2 - 

— 

-2 - 



°3 b_ 2 
ft* 3 

-G- -«= 



0.5 1 2 

ENERGY (KEV) 

Fig. 3. — Phase-resolved spectra from the summed XMM-Newton 
observations of PSR J0821— 4300. The two spectra correspond to 
two phase regions, 0.4 cycles wide, one centered on the pulse peak 
found in the softer energy band ("(pSOFT") and the other centered 
on the peak found in the harder energy band ("<f>HARD")- I n 
each case, a double blackbody model plus Gaussian emission line is 
fitted with the parameters given in Table[2] the column density and 
temperatures are linked. The cross-over energy, Across — 1.16 keV, 
of the two spectra, is shown. Residuals of the data from the model 
for the two spectra are shown in the lower panels; the soft phase 
data require a significant line component at Ey lnc = 0.79 keV. 



P\J\ — T/t c is identical to its current period P = 112 ms 
to < 1% accuracy. Thus, PSR J0821-4300 becomes 
the third CCO whose initial spin period is known to be 
large, and whose surface dipole magnetic held strength, 
B s < 9.8 x 10 11 G, is smaller than those of the canonical 
young radio pulsars. In comparison, the upper limits on 
B s fr om timing of PSR J1210- 5226 (IGotthelf fc Halpernl 
12007ft and PSR J1852+0040 ([Halpern et all 12007ft are 
< 3.3 x 10 11 G and < 1.5 x 10 11 G, respectively. 

PSR J0821-4300 is virtually a twin of the 105 ms 
pulsar PSR J1852+0040 in Kes 79 which, however, is 
~ 80% pulsed. In contrast, the weak (f p « 11%) pulse 
of PSR J0821— 4300 is more similar to the light curve of 
the 424 ms pulsar PSR J1210-5226. The latter is most 
notable for the prominent cyclotron absorption lines in 
its soft X-ray spectrum. For PSR J0821-4300, the likely 
detection of a weak emission line is reminiscent of the en- 
ergy dependence of the pulse phase in P SR J1210— 5226 
(jPavlov et al.H2002UDe Luca et al.ll2004h . which suggests 
that electron cyclotron resonance processes may be af- 
fecting its pulse shape. The surface B-field strength of 
PSR J1210-5226 is estimated to be w 8 x 10 10 G from its 
absorption features (jBignami et alJ l2003). in particular, 
the cyclotron fundamental at 0.7 keV, and the magnetic 
field of PSR J0821-4300 may be similar. 

The physics of the CCO class, evidently characterized 
by weak magnetic field and slow, original spin, may lie in 
the turbulent dynamo that generates the magnetic field. 
The strength of the dynamo depends on the rotation rate 
of the proto-neutron star ([Thompson fc Duncan! 119931 : 
iBonanno et al.l [2006ft . so the magnetic field strength is 
inversely correlated with the initial period. 



3.2. Spin- Kick Mechanisms 

Chandra proper motion measurements of RX 
J0822— 4300 reveal an extraordinar ily large trans- 
verse velocity (1570 ± 240 d 2 km s" 1 . I Winkler fc Petri 
[20071 1122 ± 360 9 km s" 1 . iHui fe Beckerl 12006b). 
As discussed by 1 Winkler fc Petrel (|2007ft . this high 
velocity restricts the possible models for NS k icks, 
of w hich there are three b road categories (see lLai I 
120011 and IWang et al.l l2006ft . The discovery of the 
spin parameters of PSR J0821— 4300 now allows us to 
eliminate two of these models: 1) the magnetic-induced 
asymmetric neutrino emission m odel requires a magnetic 
field of > 10 15 G ([Lail 12001ft and can only produce 
velocities < 2 50 km s -1 , and 2) the elec t romagnetic 
rocket effect ([Harrison fc Tademarul 119751: iLail 12001ft 
needs an initial spin period < 1 ms to achieve the 
velocity of PSR J0821-4300. Therefore, the third 
model, th at of hydrodynamic instability during the 
explosion ([Burrows et alj 120071 ). emerges as the most 
likely explanation for large NS kicks. This is consistent 
with a highly asymmetric explosion, as indicated for 
Puppis A, wherein the momentum of the oxygen knots 
match the momen tum of RX J0822— 430 0, but in the 
opposite direction (|Winkler fc Petrel 12007ft . 

Interestingly, the hydrodynamic mechanism is the one 
mechanism of the three candidate theories that does not 
lead naturally to alignment between the spin axis and 
the direction of proper motion that is observed in the 
slowly moving Crab (v t sa 120 km s _1 ) and Vela (v t « 60 
km s _1 ) pulsars, a nd is statistically ev i dent in larger 
samp les of pulsars (| Johnston et al.|[2005t iNg fc Romanil 
12007ft . If the spin of PSR J0821-4300 were aligned with 
its velocity, which is nearly in the plane of the sky, then 
we would have expected to see a large rotational modu- 
lation from any surface hot spot. On the contrary, if the 
slow initial spin period of PSR J0821— 4300 is longer than 
the duration of the kick mechanism, the impulse of the 
asymmetric kick could have been expended in one direc- 
tion rather than being averaged out by rotation. Thus, 
the long spin period of PSR J0821— 4300 allows its spin 
axis and velocity vector to be misaligned. This is consis- 
tent with our model for the weak modulation, where the 
spin axis is close to the line-of-sight, nearly perpendicu- 
lar to the velocity vector, as is allowed for a large kick 
from the hydrodynamic mechanism. 

4. CONCLUSIONS AND FUTURE WORK 

The discovery o f PSR J0821-4300 advance s the argu- 
ment, reviewed by IGotthelf fc Halpernl (|2008f ). that slow 
original spin and weak magnetic field constitute the phys- 
ical basis of the CCO class. There may even be a causal 
connection between slow natal spin and weak magnetic 
field through a turbulent dynamo that creates the mag- 
netic field. The next step is a dedicated timing study 
of PSR J0821— 4300 to determine whether it is spinning 
down steadily, and to measure its period derivative and 
dipole magnetic held. An alter native possibility that it is 
accreting from a fallback disk (Hal pern et al.ll2~007ft can 
also be tested by detailed timing. Similar studies are 
underway for the other CCO pulsars. 

We predict that most if not all of the remaining CCOs 
will turn out to be weakly magnetized pulsars with P ;> 
0.1 s. These may include the CCO in Cassiopeia A, and 
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possibly an as-yet unseen NS in SN 1987A. In the absence 
of spin parameters, is difficult to discriminate between a 
low magnetic field CCO and a quiescent magnetar on 
the basis of an X-ray spectrum alone, as these classes 
have similar spectral properties. Before the CCO pulsars 
were known, both Cas A and SN 1987A were hypothe- 
sized to host magnetars. However, the r ecent retraction 
(|Kim et all l200cfr iDwek k. ArendtJ I2008D of claimed ev- 
idence from the Spitzer Spa ce Telescope for an historic 
SGR-like outburst in Cas A (Kr ause et al.|[2005l ) leaves a 
weakly magnetized NS as the most compelling model for 
that youngest of known CCOs. In the case of Cas A, we 



must allow for the possibility that an energy-dependent 
phase shift similar to PSR J0821-4300 makes it difficult 
to detect pulsations. We are searching for similar effects 
in all CCOs. 



This investigation is based on observations obtained 
with XMM-Newton, an ESA science mission with instru- 
ments and contributions directly funded by ESA Mem- 
ber States and NASA. This work is was made possible 
by NASA XMM grant NNX08AX67X, which provided 
critical hardware support. 
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